Coupled cavity QED for coherent control of photon transmission: 
Green function approach for hybrid systems with two-level doping 
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This paper theoretically studies the coherent control of photon transmission along the coupled 
resonator optical waveguide (CROW) by doping artificial atoms in hybrid structures. We provide 
the several approaches correspondingly based on the mean field method and spin wave theory et 
al . In the present paper we adopt the two-time Green function approach to study the coherent 
transmission photon in a CROW with homogeneous couplings, each cavity of which is doped by a 
two-level artificial atom. We calculate the two-time correlation function for photon in the weak- 
coupling case. Its poles predict the exact dispersion relation, which results in the group velocity 
coherently controlled by the collective excitation of the doping atoms. We emphasize the role of the 
population inversion of doping atoms induced by some polarization mechanism. 
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I. INTRODUCTION 



Recently, all-optical and on-chip setups with the cou- 
pled resonator optical waveguide (CROW) have been im- 
plemented experimentally to show the coherent transmis- 
sion of "slow photons" [l], [2| , which is similar to the elec- 
tromagnetical induced transparency (EIT) effect [H, H[ 
occurring in the atomic ensemble medium. This success- 
ful experiment implies the coherent couplings of the sin- 
gle mode cavities [5j , which result in the photonic-crystal- 
like system with photonic band structure. For practical 
applications such as CROW setup they can be utilized to 
stop or store the light pulses propagation and then lead 
to a quantum device based on these many-body effects, 
which can also be regarded as a tunable quantum simu- 
lator for the tight binding fermion system in condensed 
matter physics. 

On the other hand, recently it was discovered [H, 0] 
that, when such an array of coupled cavities is doped 
with two-level atoms, the photon-blockaded phenomenon 
can emerge and achieve a Mott insulator state of po- 
laritons that are many-body dressed states of doped 
atoms coupled to quantized modes of optical field in the 
CROW. Most interestingly, such a hybrid system with 
a two-dimensional array of coupled optical cavities in 
the photon-blockaded regime will undergo a quantum 
phase transition from characteristic Mott insulator (ex- 
citations localized on each site) to superfiuid (excitations 
delocalized across the lattice) A similar coplanar 

hybrid structure based on superconducting circuit has 
been proposed for the coherent control of microwave- 
photons propagating in a coupled transmission line res- 
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onator (CTLR) waveguide. Here, each cavity is coupled 
to a tunable charge qubit Q . While the CTLR forms an 
artificial photonic crystal with an engineered band struc- 
ture, the charge qubits collectively behave as spin waves 
in the low-excitation limit, and these charge qubits mod- 
ify the photonic band with energy gaps to slow or even 
stop the microwave propagation in this CTLR waveguide. 
The conceptual exploration here suggests an electromag- 
netically controlled quantum device based on the on-chip 
circuit for the coherent manipulation of photons, such as 
the dynamic appearances of the laser-like output from 
CTLR waveguide where the atoms are pumped for some 
population inversion. 

These progressing investigations motivate us to further 
develop the general theoretical approach for cavity quan- 
tum electrodynamics (QED) with coupled resonators for 
coherent manipulations of photon transmission in an ar- 
tificial photonic band structure, which can be controlled 
through some new mechanisms. In Ref.[l6| we provide an 
approach based on the mean filed method and spin wave 
theory et al for different hybrid structures, which con- 
sist of the coupled cavity arrays with homogeneous (or 
inhomogeneous) couplings and various multi-level-atom 
doping. 

The present paper adopts the two-time Green function 
approach to study the coherent transmission of photons 
in a CROW with homogeneous couplings, each cavity is 
doped to a two-level artificial atom. Mathematically the 
hybrid system has the same model as that for CTLR 
waveguide connected to charge qubits Q, but the Green 
function can work well for the system which does not sat- 
isfy the low-excitation limit, in which we can even obtain 
exact solution [lfjj. We calculate the two-time retarded 
Green function for photons in the weak-coupling case. 
Its poles predict the exact dispersion relation, according 
to which the group velocity can be coherently controlled 
by the collective excitation of the doping atoms. We em- 



phasize the role of the population inversion of the toti 
doping atoms, which is induced by some polarization c 
pump mechanism. The dispersion relation exhibits son: 
exotic features such as the compressed photonic banc 
width. 

The paper is organized as followed. In section II w 
describe our setup of the photonic band device CRO\ 
interacting with doping atoms. Applying the retarde 
two-time Green function theory, in section III, we ca 
culate the eigenfrequencies of the hybrid photon-atoi 
system to characterize the coherent features of photo 
transmission. In section IV we study how the bandwidt 
and group velocity of photon transmission can be ac 
justed by controlling the doped atoms. Then we conside- 
the effects of damping in both the local mode of cavit 
and doped atoms. The stable atomic collective excit; 
tions can result in the coherent output of slow photor 
with some laser-like properties. In section V, for the phi 
nomenon of slow photons we study the effective suscept 
bility of light propagation in the CROW interacting wit 
doping atoms. In the appendices we give some necessar 
details for the Green function calculations and analy2 
the quasi-spin wave structure represented by the Green 
functions for photons and atoms that we have obtained. 



II. MODEL FOR HYBRID STRUCTURE WITH 
PHOTONIC BANDS 

We consider a hybrid structure (illustrated in Fig.l 
a )- the coupled cavity array with doping artificial 
atoms. Here, N optical cavities with homogeneous and 
nearest-neighbor couplings form a one-dimensional peri- 
odic structure, which is similar to the fermion system on 
tight binding lattice. In practice, there are two ways to 
implement such CROW. 1. With photonic crystals, the 
coupled cavities are built through regularly breaking the 
periodicity of photonic crystal. In the photonic bandgap 
materials, the cavities are defined by an array (super- 
lattice) of periodic defects in the periodic modulation. 
The inter-cavity hopping of photons is due to overlap 
between two cavity mode functions. 2. In an clcctro- 
magnetically controlled quantum device based on super- 
conducting circuit 3, the CROW is realized by the su- 
perconducting waveguide with coupled transmission line 
resonators, while the doping systems are implemented by 
the biased Cooper pair boxes. 

In Fig.l (b, c), to implement a Rabi transition with 
controllable coupling between the excited state |e a ) and 
ground state \g a ) of the doping two-level system with 
level spacing to a, the stimulated Ramann mechanism is 
usually used for a three-level system where the classical 
controlling light is resonant between the auxiliary level 
\f a ) and the excited level \e a ). 

Actually, for the sake of conceptual simplicity, here we 
assume that atom in each cavity only has three energy 
levels, two metastable lower states \g a ), \e a ) and an aux- 
iliary state |/ Q )- The transition \f a ) — > \g a ) is coupled to 




FIG. 1: (Color online) Configuration of controlled light prop- 
agation in a coupled resonator optical waveguide (CROW) by 
doping a two- level system (a). To implement a controllable 
Rabi transition between the excited and ground states of the 
effective two- level system (b), the stimulated Ramann mech- 
anism is used for a three-level system (c) where the classical 
controlling light is resonant between the auxiliary level and 
the excited level. 

a quantized radiation mode with Rabi-frcquency Q, the 
frequency uic and the creation (annihilation) operator a\ 
(a a ) is in ath cavity, while the transitions |/ Q ) — > \e a ) are 
driven by a classical controlling field with Rabi-frequency 
f2 c . Moreover, we also assume that the detuning A be- 
tween \g a ) and \f a ) with respect to the quantized light 
is the same as that between \g a ) and \e a ) with respect 
to the classical light. Due to the stimulated Ramann 
effect for large detuning A, the effective coupling can 
be obtained as g — &&l* c / A. In this sense, the effective 
coupling strength g can be well controlled by classical 
Rabi-frequency fl c and the detuning A. 

To describe the collective excitations of the doping 
atoms, we use the quasi-spin operators 

°l = |e Q )(e Q | - \g a ){g a \, (i) 
°t = l e a)(5a|, cr" = |.g Q )(e Q | 

to express the Hamiltonian H — Ha + Hac + Hq of the 
hybrid system. Here, 

JV-l 

Ha = £ f < (2) 

is the free Hamiltonian of the doping atoms and the in- 
teraction between the local atoms and the corresponding 
cavity model is of Jaynes-Cummings type 

JV-l 

Hac = 5 X! Sq(T ^ + h - c ~ ( 3 ) 

It shows a dynamic process that photons are absorbed 
when atoms transit from ground state to excited state 
while the photons are emitted when the atoms transit 



3 



from exited state to ground state. The CROW is de- 
scribed by the Hamiltonian 



N-l 



N-l 



Ha = > Loca^cia + J ^ a) a a a+ i + h.c. 



=o 



=0 



(4) 



where J denotes the inter-cavity coupling. The second 
term of He presents the tunneling of photons from the 
ath cavity to the a + 1th one. We notice that the 
model we adopt above has been used to demonstrate the 
photon-blockaded effect most recently [f| 0], and the las- 
ing behavior of the output in line resonators (CTLRs) 
by connecting each cavity to a tunable charge qubit in 
circuit QED0. 

To consider the physical significance implied by H = 
Ha + Hac + He, we perform Fourier transformations 



E 

a=0 
N-l 



z \ " °a 



Q=0 



(5) 



for k = 2im/(eN),n = 0,1, ...,7V- 1, with the peri- 
odic boundary condition for the quasi-spin operators , 
er^T = (a^y an d ° z ■ The above Fourier transformation 
describes the collective excitations of the spatially dis- 
tributed doping atoms as the quasi-spin wave via the 
collective operators a£ and a z [10| . This is because 



N-l 



°t\G) = Y.^ Ma \ E «)' 



a=0 



(6) 



represents a spin wave where \G) — \gogig2---9N-i) 
means all atoms are prepared in a ground state, while 

\ E a) = \gogi—g a -ie a g a +i—gN-i), 

means a single particle excitation in the site a. 

In order to study the collective excitation described by 
<xl and cr z , we consider the corresponding commutation 
relations 

W^k] =<7kk>, [o-fc^fc] =(jZ i 



where 



r * -i 2 - 



i N ~ 1 

z 1 \ ~* -i(k'-k)e.az 



(7) 



(8) 



means that <Tq = <r z , cr^ and er^ = can not 

generate a SU(2) subalgebra except for the case of k = 0. 
Thus and Oy can not be regarded as a collective 
angular momentum for finite N. 



Applying a discrete Fourier transformation in the k- 
space representation dfe 
tonian He, we have 



E Q e lMa a a /VN to the Hamil- 



N-l 



H = Clkajak + NuiA-n 

k=Q 

N-l 

+g E ( &k<J k + h - c -)- 



(9) 



fe=0 



The photonic band structure is characterized by the dis- 
persion relation 



Ofe = Luc + 2Jcos(k£). 



(10) 



In principle the above hybrid model can not be solved 
exactly, but we have analytically studied its realization 
based on superconducting circuit 0] when few atoms are 
populated in their excited state - the low-excitation limit. 
In this case this model will be reduced to an exactly solv- 
able coupling boson model as well as that of all-optical 
setup for stopping the light propagating in a CROW in 
Ref. @ . The crucial issue of this observation is to use the 
collective operators [10| and a~£ as bosonic spin wave 
operators in the large N limit with low excitations, since 
the usual bosonic commutation relation [c^r,cr^] = Skk 1 
can be approached as N — > oo. However, the low ex- 
citation requires a t' J k) ^ which limits the ex- 
ploitation for the general cases. So we need to develop a 
new technique to deal with the general cases. 



III. TWO-TIME RETARDED GREEN 
FUNCTION APPROACH FOR PHOTON 
TRANSPORT IN COUPLED CAVITY ARRAY 

For quantum many-body problem, the quantum or 
thermal fluctuations near thermal equilibrium may be 
characterized by time correlation functions of the type 
(A(t)B(t')), or by the Fourier transformations of these 
correlation functions, which give the correlation fluctu- 
ation spectrum. In Heisenberg picture, the time corre- 
lation function (A(t)B(t')) of two observable A and B 
depends only on the time interval t — t' by the invariant 
of time translation. The propagation of photons in our 
hybrid system can be obtained by solving Green function 
equation for photons. 

We consider the linear response of the Green func- 
tion with respect to an effective applied driving force, 
the coupling with doping atoms. We define the two-time 
retarded Green function G% (i, t') = {(A(t); B(t'))) 03 
as 

B(t'))) = -i6{t - ([A(t), B(t')}) , (11) 

where 9(t) = 1 for t > 0, and 6(t) = for t < 0. 

To obtain the equation of motion for Green function, 
we first list the equations of motion for the creation and 
annihilation operators of photons and atoms by using 



the Hamiltonian defined in Eq. ([9]). Then we have the 
equations of motion for G\ B (t, if) 



dt 



((A(t);B(t')))=-i6(t-t')([A,B}) 



-i(([A,H];B(l/))), (12a) 



^((A(t);B(t')))=iS(t-t')([A,B]) 
-i((A(t);[B,H})). 



(12b) 



F(k,k) 
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After applying the Fourier transformation from the time- 
domain to the frequency-domain, the retarded Green 
function is represented as 

({A\B)) u+iE = Jdte^+^-t'HiA^Bit'))), (13) 

where e = +0 is a positive infinitesimal. The equations 
of motion for the Green functions can be evaluated in 
frequency representation 



oj((A\B)) u = ([A, B]) + {{[A^HWB))^ (14a) 
u,{{A\B)) u = {[A,B])-((A\[B,B\))»- (^b) 

We notice that in the linear response theory, 
determines the basic spectrum structure of the hybrid 
system through the poles of ((AlB)}^ - the physical eigen- 
frequencies. 

With the notations above, we write down the equa- 
tion of the Green functions {{a k \a\)) u , ((a^, \<t^)) u , 

({akWk>))"> and ((^fc Thc Photon 

correlation ((dfe|aj,))tj can be obtained by using the com- 
mutation relation between a k and H as 



(w - n k ){{a k \a\)) w = 1 + 3«^l4»- (15) 

We can also calculate the Green function of the many- 
atom correlation ((tr^T, |cr^")) u , which satisfies 



k' 



(16) 



JV-1 



.^ e -^'- fe )^ K)+WA((f7 - k + )K 



To cut off the Green function hierarchy, we make a 
mean field avvroximation^\$i, [l4| 



(17) 



where the factor (cr„) represents the large atomic pop- 
ulation inversion in the initial state and the light-atom 
interaction hardly changes this population. The system 
of equations of Green functions has an approximately 



FIG. 2: The Feynman diagram for the effective photon trans- 
mission through the CROW in a propagating mode. The 
photon propagator contains the free part (the bare photon 
propagator is denoted by a single wiggly line) and the second 
perturbation part (one wiggly line plus a box). The box in- 
cludes the shade circle for the self energy of photon and the 
bare atom propagator( the double arrowed line). 



closed form with three simplified equations (For the de- 
tails please see the appendix A), 



, t 1 F(k,k) 

10 — ilk LO — Ilk 



(<T Z kk ,)[l + F(k,k)) 



(K"K + )>- = - 



E 

k"^k,=£k' 

V) 



Ik' M 

{cjt„ k ,)F{k",k) 



E 



fk< M ' 

<r z k , k )F(k',k) 



(18) 



fk ^ k^k fk ^ 



where we have defined 



F(fc,fc') = 52<K7K+,>> ^ 



LU — Ql 



(19) 



f k (oj) = OJ - UJ A + 



9 2 {° z ) | 
uj - Qk ' 



N-l 



\ a kk> 



To consider the basic processes of the photon distribu- 
tion in k-space, we draw the Feynman diagram (Fig. [2]) to 
interpret the above equation of ((afc|a|.)) w in terms of the 
basic processes. Here, the photon propagator l/(cu — fl^) 
(denoted by a wiggly line) appears twice. In the second 
term of the right hand side of first equation in Eq. (|18p , 
it is modified by an interaction with atomic flips char- 
acterized by bare atom propagator ((c^T |cr^"/)) u , which is 
denoted by a double line. This Feynman diagram de- 
scribes a second order process of the interaction between 
the localized modes of the optical field and the doping 
atoms. 

We consider the weak-coupling case. On the right side 
of the equation of ((a^ \<J k })ui in Eq. (fi"8|) there are two 
terms, one is about the zero order of coupling constant 



5 



g and the other is about the second or higher order of g. 
Here we solve the 
order term of g: 



Here we solve the equation of ((a, \<J^)) UJ with the lowest 



fk M 



So we obtain the lowest order solution for 



({ak\al)) u = 



U) — UJa 



and 



(K"k + »o 



(W - - fife) + 5 2 ((7*)' 

-(a z )(^-»fe) 

(w - n fe )(w - wa) + .g 2 (o- 2 ) ' 



(20) 



(21) 



(22) 



There exist two poles lu — and lu — tu k 

w^ 5 = flu ± e fc , (23) 
which are the dispersion relations. Here 
1 



(24) 



We analyze the properties of the retarded Green func- 
tions of ((afe|a fe )) w+ie and ((cr^|(T^)) a , + i e by decomposing 
them into two branches of wave, respectively. The prop- 
agating photons have the form as 

({a k \al)) u+is = A k G+(k,u) + B k G-(k,u), (25) 

where the free Green functions are denoted by 

1 



G ± {k,uj) 



(±) , • ' 

IjJ - 0J k + IE 



(26) 



and the atomic part has the form as 

«0*>£» W -H« = -(<J z )[B k G+(k,uj) + A k G-(k,uj)}, 
where 



(+) (-) 

A k = — — , £>fc — 



(-)■ 



W fc _ U k 



(27) 



are the amplitudes of the two wave branches. 

Transforming the Green functions above back to the 
real time representation, one can observe that the pho- 
tons propagate with two frequencies 0J k ^- If we regard 
the transmission of the localized photons as propagating 
wave, the two wave branches or two partial waves with 
iJu^ have the amplitudes A k and B k , respectively. For 
(( a k \ a t ))f> ^ nas been observed that [1] the total collec- 
tion of the identical two-level atoms can be regarded as 
an ensemble of N spins and thus its collective excitation 
can be described as spin waves, which are characterized 
by (( a k \ lJ k ))u- It is shown that the spin wave has two 



eigenfrequencies , but with the twist amplitudes B k 
and A k , respectively. (The detailed analysis is given in 
the appendix B.) 

For the photons in the CROW from localized modes to 
propagating modes, we can now visualize the two wave 
branches as quasi-particle excitations by considering the 
existence of isolated poles of ((6fc|a£))w Suppose that 
there are such poles to k — > u> k — vy± on the complex 
plane. They correspond to the life l/7± of the quasi- 
particle excitations characterizing the two branches of 
propagating wave in the coupled cavity array. By phe- 
nomenologically adding imaginary parts — i~/ c and —ija 
to the cavity eigenfrequency uic and the atomic level 
spacing lua respectively, 7± can be explicitly expressed 
obviously in terms of —i"f c = — and —ija = — *7 (the 
details in the next section). This means that the decay of 
transporting photons is just induced by the cavity decay 
and the atom natural line width. 



IV. COHERENT TRANSMISSION OF 
PHOTONS WITH SLOWED GROUP VELOCITY 



From the dispersion relation ([24]) . it can be observed 
that the population inversion (<j z ) can directly affect the 
basic features of coherent transmission of photons in the 
CROW. To enhance this influence, we put more doping 
atoms in a cavity. Suppose that every cavity is doped 
by n identical atoms without interaction among them- 
selves. Then the parts of the Hamiltonian in section II 
concerning atoms become 



N-l 



N-l 



a=0 



a=0 



where we have introduced the collective spin S a = 

S™=i <TqJ - I n this sense the above frequencies w£ can 
be modified by replacing (a z ) with 



JV-l 



(28) 



Q=0 / 



Obviously, the average of the total spin is bounded as 
-n < (S z ) < n. 

Before discussing the group velocity of photons, we 
investigate the change of bandwidth of this photonic- 
crystal-like system. Because the group velocity v k g can be 
calculated according to v k g = duj/dk, which concerns the 
range of lu, the change of bandwidth plays an important 
role. Without the doped atoms the spectrum of photons 
should have only one band, and the central line should 
be at loc- However when atoms are doped, the spectrum 
splits into two bands with eigenfrequencies u k ■ Then 
the central lines shift to £Id ± e k=n /2^- Without popu- 
lation inversion, i.e.,(5 z ) < 0, the two bands have the 



same width W — 



,(±) 



|fc=o 



,.(±)i 

w fe \k=n/2e 



< 2 J, which 
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is calculated as 



W = \F ± ((S$),0)-F ± ((S* ),-)\, 



(29) 



where 



F± (x,k) = J^{8 + J cos(H)] 2 - g 2 x, (30) 

for 8 = ujc~uja- This means that the bandwidth becomes 
narrower when cavities are coupled to more atoms. 

Next we consider the group velocity of photon propa- 
gation 



w (±) fc = J£ S m(k£) 



1 ± 



5 + 2Jcos(k£) 



(31) 



2F±«S*),fc) 

for various cases. At k = tt/2£, the group velocities of 



W £/L read 




0.02 0.04 0.06 0.08 0.1 
g/J (a) 



0.02 0.04 0.06 0.08 0.1 
g/J (b) 




0.02 0.04 0.06 

g/J (c) 



0.02 0.04 0.06 0.08 0.1 
g/J (d) 



k 



(32) 



and the amplitudes of the photon propagator can be cal- 
culated as 



.4 



8 + k 
2k 



1 — 



2k 



(33) 



FIG. 3: There are ten curves with J > S > 0.1J in the 
fig(a)-(d), respectively. The group velocities of Vg, Vg~^ and 
the amplitudes Ah, are functions of coupling strength g 
at k — n/(2£). In the fig(a)-(d) the lattice constant t = 1, 
| < S z > | = 10. In fig(a) and fig(b) the upper curves are 
for S — J, while in fig(c) and fig(d) the upper curves are for 
5 = 0.1J. 



respectively, where x = 8 2 — 4g 2 (S z ). We now con- 
sider the case with most atoms in the ground state, i.e., 
(S z ) < 0. When 8 > 2gy/\(S z )\, the amplitude at the 



band center A 



■jr/21 



1, B 



TV /2l 



((a 7I /2£\al /2e )) 



0, and then 
1 



uj-\-ie 



\j2l 



(34) 



IE 



In other words, in this limit the photon modified by 



(+) 

/2V 



Cc 



the atoms tends to have an eigenfrequency or 
respondingly the group velocity reaches its maximum 
Vg" ~ 2 J£, and the quasi-spin wave for the atomic 

excitations is characterized by the Green function 

\{S*)\ 



i( (7 Tr/2l\ <J Tr/2£))u+i£ 



\/2t 
(") 



(35) 



IE 



which has a distinct eigenfrequency wV,^' With g — > 0, 
wi^o» and a;i/^ approach and o;,a, respectively. On 



k/21 



j -k/21 



the contrary, when 8 <C — 2<7y|(iS*)|, we make a similar 
argument: as A v /2i ~^ 0, 13^/21 - > 1, the photons and 
atomic spin wave propagate with eigenfrequencies ^\\n 

and 



and ^>^\v res P ec tively. By letting g 

w i/2£ can r,e rev i vecl as W A ancl wp, respectively. The 
group velocity of photons can also reach its maximum 



0, u 



(+) 



-k/21 



(-)ir/2e 



2J£ (The detailed analysis is given in the 



appendix B). These observations are different from the 
results obtained in the simple cavity - cavity coupling 
system without atom doping in Ref.[5[ . 

Analyzing the features of eigenfrequencies for 
photon and atom parts in the case of weak coupling, we 

have observed that have different preferences to ap- 
proach frequencies of pure photons or bare atoms. It is 
concluded from this observation that , if 6 > 0, photons 
prefer wj^ while the atomic spin wave prefers w[. ■* ; if 
8 < 0, the conclusion is just on the contrary. We illus- 
trate these analysis in Fig. [4] 
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FIG. 4: In fig (a) and (b), a/^^j (solid) and u^/L e \ (dashed) 
are functions of coupling strength g and j < S z > J = 10. In 
fig (a) u>c = 2J, loa = 1.5J and 5 — 0.5J, while in fig (b) 
luc = 1-5J, loa = 2J and 8 — — 0.5J. 
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Next we consider how a coherent pump induces pop- 
ulation inversion to result in a laser-like output for the 
CROW. In the discussions above, we have considered an 
ideal case in which the quantum dissipation and dephas- 
ing due to the influence of the environment are neglected. 
Meanwhile, we have not concerned the role of (S z ), the 
average value of total atoms population. Here, we also 
assumed that (S z ) is not tunable . But for an open sys- 
tem, (S z ) becomes a time-dependent parameter. We can 
tune the population of atoms to change the properties 
of photon transmission. When the population inversion 
takes place, we expect that laser-like output emerges. 

To see more details, let us consider a realistic case that 
the cavity damp has the same rate K and the atoms have 
the decay rate 7 due to spontaneous radiation. Then 
luc — * LL>c' — i K an d to a — * u A — *7, so the eigenfrequencies 
of photonic band become 

u; k ' =il D ± to±, (36) 

where we define A = 7 — k and 

Ob = \^m k -Lo A )-iKY-^{S z ). (37) 

For the sake of simplicity, we consider a special case 
that is £l k « u>a, and 7 k. Then the eigenfrequencies 
read 

4^ =n D ± l -{f((S z )hTjTK], (38) 

where f(x) = yl + 4xg 2 /^ 2 . If the imaginary part A of 
an eigenfrequency (e.g., w£ ) ^ s positive, the laser-like 
output will appear since a component Af./(u> — w£ ) of 
{{ak\a\.)) u has a real time correspondence 

A(t) = f Akdu) ~ -ie(t)e- iQot + Xt , (39) 

where A k is a slowly varying amplitude in the fc-space. 
Actually, when most of atoms stay at the ground state, 
i.e. (S z ) < 0, it is impossible for and to have 
a positive imaginary part; when the population of most 
atoms inverts, i.e. (S z ) > 0, the imaginary part of lo^ 
may be positive. 

Furthermore there is a threshold value of (S z ) to satisfy 
the condition [/((#*))- 1] 7 > k. When (S z ) < 7 2 /(4 5 2 ), 
we explicitly obtain the threshold value of (S Z ) T as 

(Sir = (40) 

Above this threshold value, the eigenfrequency 

J k + ^n D + l (£(s z )-^) (4i) 

has a positive imaginary part, which results in a laser- 
like output in the CROW. It is very interesting that the 



(S z )t is very similar to the threshold value of popula- 
tion inversion in the generic laser theory. We also notice 

that, in weak-coupling limit, u>i ^ can not be a robust 
frequency of laser-like output since it damps fast with 
rate [/ + 1] 7 + «. 



V. SUSCEPTIBILITY ANALYSIS FOR LIGHT 
PROPAGATION IN THE DOPED CROW 

The above analysis displays a possibility to implement 
a slow light propagation in the doped CROW, but the 
calculation of the group velocity from the dispersion re- 
lation shows that, only for certain wave vector k, can the 
group velocity be reduced down. Thus for the propaga- 
tion of a wave packet or a light pulse, we still need some 
details for the absorption and dispersion of light in the 
doped CROW. We use the dynamic algebraic method de- 
veloped for the atomic ensemble based on quantum mem- 
ory with EIT [1 01 ] . The original method was proposed for 
the conventional EIT system, which consists of a vapor 
cell with three-level A— atoms near resonantly coupled to 
the controlling and quantized probe light. Our dynamic 
symmetry analysis is based on the hidden dynamic sym- 
metry described by the semi-direct product of quasi-spin 
SU(2) and the boson algebra of the excitations. This 
method allows us to build a dynamic equation describing 
the propagation of the probe light in this atomic ensem- 
ble with atomic collective excitations fill ] . 

Now we apply this algebraic method to calculate the 
susceptibility of light for the group velocity of photonic 
wave packet propagating along the doped CROW. Then 
we investigate how the susceptibility depends on the var- 
ious control parameters. 

We simplify our model by using the collective oper- 
ators b a = S~/y/n and &t — S+/y/n, which represent 
the quasi-spin wave in the low-excitation limit that only 
few atoms populated in their excited state. In this case 
we can check that the spin wave is bosonic excitation 
since the boson commutation relation [6 a ,6^,] = S a a i 
is satisfied in the low-excitation limit. Then the total 
Hamiltonian for the hybrid system with many-atom dop- 
ing becomes 

a a 

+ aj A J2 b ^ + ^29Vn(alb a + bla a ) . (42) 

a a 

Its fc-space representation H = J^k is a simple sum 
of the fc-component 

H k = Vl k a\a k + uj A blb k + g^/n ^a^.6 fe + h.cj . (43) 

Here, we have used the Fourier transformation b k = 
J2 a exp(ikea)b a /VN. 

For each mode k, we can write down the Heisenberg 
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equations of operators a k and b k . 

id t a k = -iKa k + Q k a k + gy/nb k , 
id t b k = -ijb k + u) A b k + g\fna k . 

Here, we have phenomenologically introduced the decay 
rates k and 7, and 7>k. In the interaction picture, we 
adopt the time-dependent transformation, 



a k = a k e inkt ,b k = b k e 



(44) 



for a k and b k to remove the fast varying parts of the 
light field and the atomic collective excitations. Then 
the above equations of motion are reduced into 

d t a k = -na k - ig^/nb k , 

d t b k = -jb k - i{u) A - tt k )b k - igy/na k . 

In general the steady state solution of the equations 
above determines the susceptibility of photon transmis- 
sion. It is noticed that the quantized light described by a k 
is the superposition of some localized modes aj . On the 
contrary, the spatially distributed photon field is charac- 
terized by a a = ^2 k exp(—ik£a)a k /\/N, which means the 
inhomogeneous polarization (P a ) depends on the spatial 
position. Correspondingly, we have the fc-space represen- 
tation of the light field 



E k (t) 



a k e 



—iQkt 



h.c. 



(45) 



In comparison with the classical expression E k (t) 
e k exp(— if^i) + h.c, it is recognized that 



a k . 



On the other hand the linear response of medium 
is described by the local polarization (P k ) = 
(p k ) exp(— iQ k t) + h.c, where the average polarization 



(pk 



Vri(b k 



(46) 



slowly varies and determined by an average value of exci- 
tation operator b k ; /i denotes the dipole moment of single 
atom, and V is the effective mode volume[l5j]. It is re- 
lated to the susceptibility \k of the /c-space by 



Xk = 



(Pk) _ Vny 
(ek) eo (fife) e Q V 



(47) 



since (P k ) = XkE k (t). 

To calculate the susceptibility in our case, we need the 
steady state solution satisfying d t b k = 0, or 

7&fe + i(— S — 2 J cos kt)b k + ig^/na k = 0, 

for 7»k. Here, S = loc — is the detuning between 
photons and atoms. In the steady state approach, we can 
take the expectation value for the above equation 

ig\/n (a k ) = — (7 — iS — 2 J cos k£)(b k ) . 



Since the dipole approximation g = —\i ^luc / C^Ve®), the 
linear susceptibility 

Xk = Xik + iX2k 

2ig 2 n 



The real part 



Xik 



uj c [i - i(S+ 2Jcos ki)\ 



-((5 + 2Jcos kl)2g 2 n 



uc\l 2 + (S + 2J cos kl) 2 Y 
and the imaginary part 

2g 2 n r y 



X2k 



oj c [l 2 + (5 + 2J cos Uf 



(48) 



(49) 



are related to the dispersion and absorption of the light 
field in the CROW, respectively. 




3-2-10123 
6 (b) 





-3-2-10123 
6 (d) 



FIG. 5: Real part xi (solid) and imaginary part X2 (dashed) 
of the susceptibility \ vs the light detuning 5 in normalized 
units of 7 according to: (a) k — ir/2£ and J = 0.1; (b) k = n/£ 
and J = 0.1; (c) k = and J = 0.1; (d) k = and J = 0.8. 
The other parameters are given as: loc = 1 <?\A* = 1- 

Because the photons have the band structure, the 
properties of dispersion and absorption of photons vary 
with the the wave vector with momentum index k. While 
the real part xik reaches its maximum at 5 = — (7 + 
2 J cos W), the imaginary part \2k reaches its maximum 
at S = —2Jcosk£ and thus the absorption is a consid- 
erable property of the system. Here, the photons with 
different wave vectors k will have different character of 
absorption. Figure [5] (a)-(c) show the dependence of 
Xi k and %2fe on k. The maximums of the absorption 
for k — ir/2£, ir/£, 0, appear at three different values of 
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d. The reason for this phenomena is that the inter-cavity 
interaction with the coupling constant J will shift the res- 
onance point in general, but if k = n/2£, the coupling has 
no effect on for spectral structure. In view of our anal- 
ysis, the absorption directly depends on the wave vector 
k. At the same time we can imagine that the character of 
absorption can influence the group velocity of photons. It 
is obvious to see that an unavoidable loss effect appears 
for the group velocity when the atom media absorbs light 
strongly. The dispersion relation of photons is described 
by f2fc = iuc + 2 J cos ki , from which we calculate the 
group velocity as 

^'Hf |cx2J£ - (50) 

But from Eq. (j4"9")) it can be seen that the media absorp- 
tion characterized by X"2k, will be stronger when 2J£ be- 
comes smaller. In other words, the considerable absorp- 
tion corresponds to a slow group velocity. Actually, since 
the effect of the group velocity is due to a spectrum struc- 
ture of the wave vector fc, only a small range of k around 
the point corresponding to the minimum group velocity, 
avoided is the higher order dispersion. The the point of 
minimum group velocity and the point of maximum ab- 
sorption are related and somewhat close. This fact means 
that there is some unavoidable loss. 



VI. CONCLUSION 

We have studied the coherent transmission of photons 
with local modes along the CROW coupled to artificial 
two-level atoms. Under the weak-coupling limit, we use 
the stimulated Raman excitation to tune the level spac- 
ing of the effective two-level system so that the properties 
of photons in the CROW can be manipulated coherently. 
As the above results display, if we prepare the hybrid 
system as that loq 3> ^>a or ujc *C uja, the group veloc- 
ity of photons in the doped CROW will reach maximum 
under the two cases. Meanwhile the two eigenfrequen- 
cies of the hybrid system have preference that, while one 
tends to the frequency of photons, another tends to that 
of quasi-spin wave of the total atoms. By controlling the 
average population of the doping atoms in the CROW 
with decay, we predict that the laser-like output may oc- 
cur. With such an exotic photonic band structure, the 
light with different k has different properties of absorp- 
tion. 

In about control of photon transmission in CROW 
by doping artificial atoms for various hybrid structures, 
we study the case of the resonate three-level doping 
atoms by making use of the quasi-spin wave theory based 
on a mean field method. This investigation will make a 
corporate effort for the coherent transmission of photons 
in an artificial structure, where both the EIT effect and 
the band-like structure are utilized simultaneously. 

This work is supported by the NSFC with grant 
No.90403018, 10375038, 60433050, and NKBRPC 



with No.2004CB318000, and NFRPC with No. 
2001CB309310 and 2005CB724508. One (LZ) of 
authors gratefully acknowledges the support of K. C. 
Wong Education Foundation, Hong Kong. 

APPENDIX A: EQUATIONS OF MOTION FOR 
THE TWO TIME GREEN FUNCTION 

In this appendix, we provide a detailed derivation of 
the approximately closed system of the equations for two- 
time Green functions as following text 

({ak\a{)) u , ((ovlcr^))^ ({<7kWt))m 

«a*k£»«> <Kl4'»- ( A1 ) 

From the commutation relation between a k and H, the 
equation of ((ak\a\.))u) is obtained as 

u{{a k \al)) u = ([a k ,al]) + (([oA,,-H]|4»a;. 

or 

(w - n k )((a k \al)) u = 1 +ff«o£-|o£)) u . (A2) 

Since the above equation concerns the two-time Green 
function ((& k \a> k ))u)} we need its motion equation, but 
here, we first calculate the equation oi{{(Ty \cr k ))u 

or 

"((viWt »" = UA((<ri\°k))» ( A3 ) 
-^2^ e (\° a a k"W k ))«> 

The mean field approximation assumes that (cr*) can be 
factorized from ((a^a k " i.e., 

<<^M4>>" « <<><<M4>>" (A4) 

and then the above Green function hierarchy is cut off. 
Thus we get a system of Green function equations 

w«<vi<7+»„ « -|E e " (fc " fc " )£Q «)«^"i^))- 

k"a 

-1 N ^ 

e^ k - k)la {al)+u A {{° k ,\o-t)U 

or 

v ' a— 

+ gJ2e- i{k '- k " )ea (o- z a )(Ma+)U. 

k" a 
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For k — k' , the Green function ((<J k Wt)) u satisfies and 



JV-l 



a=0 

gJ2e- i{k - k ' )ea K)((a k >\a+)U. (A6) 

k' a 

We notice that the equation of {(Sfc|0fc/))w is given by 

^«a fc |<r+)) a; = (([a fc) ff]|4)) w , 

or 

(w - r! fc )((a fc | ( 7+)) w = . 9 (K7K+}}.. (A7) 

In order to derive the equation about ((o' k \a, k ,)) U j , we 
use the second kind of motion equation (|14bj) . 

"«*fcl4»" = -{fa\[&l>E\))u> 

or 

(u-n k ,)((a^\al)) u =g((a^\a+)) u . (A8) 
By defining 

N—l 



'k'kl 



TV ^ 



e-^ k - k ^ a (a z a ), 



f k (U>) = (UJ -Ll> A + 



9 2 (<? z ) ■ 



the equations about ((a k \<J k ))u and ((c fe/ |a^!"))aj can be 
finally obtained as 



/ fe H(K7| ( 7+)). = -< ( T z )-.g 2 ^K, 



(K<K + » 



(A9) 



/*' H((^K + )) w = -fe) 



(A10) 



2 V- ,z \«0-fc«kfc»W 

5 2^ (^"fc' 



k"^k^k' 



LU - fifc// 



APPENDIX B: QUASI-SPIN WAVES COUPLED 
TO TRANSFERRED PHOTONS 

In this appendix we analyze the physical meaning rep- 
resented by the Green functions for photons and atoms 
that we obtained in section III. First we explicitly rewrite 
the coefficients A k and B k in ((ofc|o£)) w and {{(T k \(J k ))w 
as 



A k = 



(Bl) 



(B2) 



where A& = fifc — w^. Let wp-Wi = 5. When k = ir/(2£), 
we obtain G P = ((a 7I / 2 e\al /2e )) UJ as 



G P = 



6+ yJS 2 -4 9 2 (<t z ) 



1 



-<5+ v^ 2 -V(cr z ) 1 

FT' 



-V (a* ) cj - - 



(B3) 



From the above equation we can see that, when (a z ) < 
0, and 5 ^> 2g-\f\{a z )\, the amplitudes at the band center 



.4 



n/2l 



1) #71-/ 



2/' 



0, 



which means 



({K/2t\% /2 ^))u 



.(+) : 



(B4) 



and the group velocity w. 



2J£. Meanwhile, if 



g — ► 0, the eigenfrequencies of photons and atoms corre- 
spondingly approximate to their original eigenfrequencies 
without coupling 



,(+) 



J ^/2l^ ^A- (B5) 

If detuning 8 <C — 2g-y/[(a^)|, the values of amplitudes 



arc in reverse 



t/2£ 



0, -Btt/2^ 



1. 



(B6) 



Thus the Green function of photon at the band center 
only has one wave 



(vV2*|a]r/2*»« 



1 



"V/2« 



(B7) 



It also can be obtained that v g « 2J£. Meanwhile, 
we can conclude that, when g — > 0, the eigenfrequen- 
cies of photon and atom are recovered correspondingly 
by another way that 



(B8) 



Next we study the Green function of the doping atoms 
GA=(Kk+))„: 



11 



with amplitudes 



and 



Bl 



,(+) 



-n* 



,.(+) ,,(-) 



n k - 



k - & 



(-) 



4 +) 



= A k 



(BIO) 



(Bll) 



which has the similar expression as those of photons. 
Thus we rewrite the atomic Green function as 



G A = -(<r z 



A u 



,(+) 



,(-)■ 



(B12) 



We also consider the situation at k = tt/(2£). First 
we assume (a z ) < and S ^> 2gy/\(a z )\, in this case, 
the value of the amplitudes approximate to one and zero 
respectively 



A 



tt/2£ 



1, B„/ 2 



2f 



0, 



(B13) 



and thus the Green function of the doping atoms becomes 

1 



\/2£ 



(BU) 



However, when g — > 0, 

,(-) , 



,(+) 



If (a z ) < and S < -25r v /|(cr 2 )|, we have 



-4, 



/2* 



0, -B,r/2£ — > 1, 



and thus 



(K/2eWi/2e))«> ~ IK) I ?+f 



w/2( 



Meanwhile when g — ► 0, the eigen-frequencies 



(B15) 



(B16) 



(B17) 

Finally we conclude that if 5 3> 2g-\/|(c7 z )|, ^ s 

the eigenfrequency of photonic part while w ^X^ is the 
eigenfrequency of atomic part; On the other hand if S <C 
—2gy/\(cr z )\, the w^7^ is the eigenfrequency of photonic 

part and w^t^ is the eigenfrequency of atomic part. 



[1] Q. F. Xu, S. Sandhu, M. L. Povinelli, J. Shakya, S. H. 

Fan, M. Lipson, Phys. Rev. Lett. 96, 123901 (2006). 
[2] R. W. Boyd, D. J. Gauthier, Nature 441, 701 (2006). 
[3] S. E. Harris, Phys. Today 50, 36 (1997). 
[4] M. D. Lukin, A. Imamoglu, Nature (London) 413, 273 

(2001). 

[5] M. F. Yanik, S. H. Fan, Phys. Rev. Lett. 92, 083901 
(2004). 

[6] D. G. A ngelakis, M. F. Santos, S. Bose, 

|quant-ph/06"06159 (2006). 
[7] M. J. Hartmann, F. G. S. L. Brandao and M. B. Plenio, 

Nature Physics 2, 849 (2006).. 
[8] D. Greentree, C. Tahan, J. H. Cole, L. C. L. Hollenberg, 

Nature Physics 2, 856 (2006). 
[9] L. Zhou, Y. B. Gao, Z. Song, C. P. Sun, 



|cond-ma t/0608577 (2006). 
[10] C. P. Sun, Y. Li, X. F. Liu, Phys. Rev. Lett. 91, 147903 
(2003). 

[11] Y. Li, C. P. Sun, Phys. Rev. A 69, 051802(R) (2004). 
[12] D. N. Zubarev, Soviet Phys Usp(English Trasl), 3, 320 
(1960). 

[13] P. A. Wolff, Phys. Rev. 120, 814 (1960). 

[14] T. Izuyama, D. J. Kim, R. Kubo, J. Phys. Soc. (Japan) 
18, 1025(1963). 

[15] M. O. Scully and M. S. Zubairy, Quantum Optics (Cam- 
bridge University Press, Cambridge, 1997). 

[16] Lan Zhou, Jing Lu, C. P. Sun, quant-ph/0611159, sub- 
mitted to Phys. Rev. A. 



